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ABSTRACT 


Perturbed  acoustic  propagation  is  considered  between  a point 
source  and  an  array  receiver  in  which  the  perturbations  are  due  to  a 
Garrett-Munk  internal  wave  field.  The  resultant  acoustic  phase  fluctua- 
tions along  the  array  are  determined.  Straight  ray  propagation  is  assumed, 
and  three  array  orientations  are  considered:  broadside  horizontal,  end- 
fire  horizontal  and  broadside  vertical.  It  is  shown  that  for  an  acoustic 
frequency  of  150  Hz,  a range  of  50  km  and  a horizontal  broadside  receiver 

separation  distance  (in  the  array)  of  1 km,  the  rms  phase  difference  between 

$ 

the  two  receivers  at  1000  m depth  in  the  N.  Pacific  is  typically  5 , and  in 

t M 

the  Arctic  it  is  typically  2r  , whereas  in  the  N.  Atlantic  at  the  same 
depth  it  is  typically  50^  . These  geographical  variations  are  due  mainly 

A 

to  variations  in  the  potential  sound  velocity  gradient.  Similar  phase 
differences  occur  for  a vertical  receiver  separation  of  50  m.  The  hori- 
zontal receiver  separation  necessary  to  reduce  the  acoustic  coherence  to 
0.5  for  the  same  operating  conditions  is  shown  to  be  very  large  in  the 
N.  Pacific  (3  400  km)  and  Arctic  (>  1000  km),  but  only  2.3  km  in  the 
N.  Atlantic. 
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INTRODUCTION 

Any  system  which  makes  use  of  phase  comparisons  in  underwater 
acoustic  signals  will  ultimately  be  limited  in  its  accuracy  by  the  random- 
izing effect  of  the  natural  state  of  motion  in  the  ocean  medium.  In  recent 
years,  this  state  of  motion  has  been  the  subject  of  much  research,  parti- 
cularly for  those  length  and  time  scales  than  can  be  described  as  internal 
waves.  One  of  the  results  of  this  research  is  that  there  now  exists  an 
analytic  model  than  can  be  used  to  predict  acoustic  properties  of  the 
ocean,  such  as  the  expected  rms  acoustic  fluctuations,  the  coherence 
lengths,  frequency  and  wavenumber  spectra,  etc.  The  internal  wave  model, 
known  as  the  Garrett-Munk  spectrum  (Desaubies,  1976a  and  b),  is  empirically 
based  and  has  been  subjected  to  some  experimental  confirmation  since  its 
proposal  (for  a general  reference  see  the  collection  of  papers  edited  by 
Briscoe  (1975)). 

If  an  acoustic  signal  is  transmitted  over  a fixed  path  between 
two  points,  the  phase  difference  between  the  received  and  transmitted 
signals  will  be  dependent  on  a number  of  variables,  for  example:  signal 
frequency,  path  length,  the  total  sound  velocity  structure,  receiver  and 
source  depth,  etc.  If  the  sound  velocity  structure  is  unvarying  in  time, 
a signal  at  a given  frequency  will  exhibit  an  unvarying  phase  difference. 

If  the  sound  velocity  structure  changes  in  time,  the  given  signal  will 
exhibit  varying  phase  differences,  and  the  characteristics  of  the  varia- 
tion will  be  intimately  linked  to  the  characteristics  of  the  velocity 
changes.  If  the  sound  velocity  changes  are  due  to  background  internal 
waves,  they  will  be  distributed  along  the  entire  acoustic  path  and  the 
resultant  phase  difference  the  receiver  will  be  an  integral  of  these 

changes  taken  along  the  path.  In  the  treatment  given  here  the  acoustic 
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phase  perturbation  is  defined  as  the  extra  amount  of  phase  difference  caused 
solely  by  the  presence  of  the  internal  waves. 

The  internal  wave  field  is  modelled  stochastically,  that  is, 
only  expected  values  are  known,  and  the  internal  wave  field  at  any  given 
instant  is  merely  one  unspecified  realization  taken  from  an  infinite  ensemble 
of  such  realizations.  The  acoustic  phase  perturbation  is  thus  a random 
process  and,  therefore,  only  expected  values  can  be  predicted. 

Two  papers  have  recently  been  published  that  deal  extensively 
with  this  acoustic  perturbation  problem.  They  provide  overlapping  treat- 
ments and  so  they  will  be  discussed  here  together.  The  first  is  by  Munk 
and  Zachariasen  (1976),  hereafter  referred  to  as  MZ,  the  second  is  by 
Desaubies  (1976a).  In  these  papers,  the  authors  have  shown  that  for  the 
two-point  fixed-geometry  single-receiver  acoustic  transmission  situation 
referred  to  above,  the  expected  rms  phase  perturbation  caused  by  internal 
waves  is  proportional  to  the  signal  frequency  and  approximately  proportional 
to  the  square  root  of  the  path  length.  The  proportionality  constants  are 
governed  by  parameters  from  the  internal  wave  model.  They  have  also  shown 
that  the  perturbations  change  very  slowly  in  time;  the  frequency  spectrum 
of  the  phase  perturbations  is  heavily  weighted  towards  the  inertial  period 
(15.7  hrs  at  50°N  and  12  hrs  at  90°N). 

For  this  single-receiver  configuration  then,  and  assuming  no 
other  perturbations  are  present,  the  received  phase  will  appear  to  be  stable 
over  short  time  periods  (say  \ hr  or  less),  although  the  value  of  the  received 
phase  will  be  different  from  the  unperturbed  value  by  a random  amount  with 
a root-mean-square  that  is  governed  by  the  internal  wave  field  parameters 
and  details  of  the  configuration. 

For  an  array  receiver  configuration  and  a single  source,  Figure  1, 


the  received  phase  at  any  single  hydrophone  will  also  appear  to  be  stable 
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Figure  1.  Schematic  representation  of  the  propagation 
conf igurat ion. 


over  short  times  but  the  difference  in  phase  between  two  separated  hydro- 
phones will  exhibit  spatial  variations  which  are  caused  by  the  spatial 
inhomogeneities  in  the  internal  wave  field  realizations.  The  phase  differ- 
ence that  exists  simultaneously  at  spatially  separated  hydrophones  will  have 
a mean-square  (in  the  ensemble  sense)  that  approaches  zero  at  very  small 
separations,  and  generally  increases  as  the  separation  increases. 

If  the  joint  probability  density  of  phases  at  the  two  hydro- 
phones is  Gaussian,  and  the  received  amp  1 itudes  are  unperturbed,  the 
coherence  between  the  two  received  acoustic  signals  will  be  exp { — %<t> 2 > , where 
<t>!  is  the  mean-square  phase  difference.  Thus,  the  acoustic  coherence  will 
be  less  than  0.5  whenever  the  rms  phase  difference  is  larger  than  67.4°. 

Of  course,  in  a real  propagation  situation,  the  acoustic  amplitudes  are  also 
expected  to  be  randomized  by  the  internal  wave  field  with  the  result  that 
the  acoustic  coherence  will  fall  below  0.5  for  smaller  values  of  than 


67.4  . Unfortunately,  the  treatment  of  acoustic  amplitude  perturbations  is 
much  more  complicated  than  that  for  acoustic  phases  (see  MZ  who  show  that 
diffraction  effects  must  be  included  for  the  former  but  not  for  the  latter) 
and,  in  the  interest  of  expediency,  the  analysis  given  in  the  following 
section  will  treat  only  phase  fluctuations. 

The  purpose  of  the  following  analysis  is  to  provide  estimates  of 
the  expected  amount  of  rms  phase di  fference  for  two  receivers  situated  in 
typical  N.  Pacific,  Arctic  and  N.  Atlantic  conditions.  Three  array  orien- 
tations will  be  considered:  endfire  horizontal,  broadside  horizontal,  and 
broadside  vertical.  The  numerical  estimates  are  obtained  by  using  straight 


ray  approximations  and  the  Garrett-Munk  spectrum. 

The  receiver  separation  necessary  to  reduce  the  acoustic  coherence 
to  0.5  (i.e.  rms  phase  difference  67.4°)  is  calculated  for  the  above 
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conditions  for  a horizontal  broadside  array  at  a range  of  50  km;  and  predic- 
tions are  also  made  (in  Appendix  A)  using  a variant  of  the  Garrett-Munk 
spectrum  due  to  Desaubies. 

Surface  and  bottom  reflections  are  specifically  excluded  from 
all  calculations. 

FLUCTUATION  MECHANISM 

There  are  two  predominant  mechanisms  by  which  internal  waves  can 
affect  acoustic  signals:  (i)  raising  and  lowering  of  surfaces  of  constant 
sound  velocity,  and  (ii)  Doppler  shifts  and  advection  due  to  fluid  particle 
flow.  It  is  shown  by  MZ  that  the  former  is  an  order  of  magnitude  more 
important  than  the  latter  so  only  the  former  will  be  treated  here.  The  wave 
equation  for  sound  propagai.  'n  in  tn>'  absence  of  advection  is 


(1) 


and  in  the  presence  of  internal  waves  (or  any  perturbing  agency) 


'(c° 


(z)  + 6c ( x , y , z , t ) 


c 2 ( z ) + 2c  ( z ) 6c 
o o 


= c2 (z ) (1  + 26c/c  } 
o o 


(2) 


In  these  equations  p is  acoustic  pressure,  c is  sound  speed  and  c^  is 
unperturbed  sound  speed. 

It  is  shown  by  MZ  and  Desaubies  (1976a)  that  the  effect  of  the 
inhomogeneities,  6c,  on  the  phase  of  p can  be  calculated  by  ray  tracing  or 
ray  timing  techniques  (at  least  for  paraxial  rays). 
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From  the  ray  model,  the  travel  time  T from  source  to  a single 
receiver  is  given  by 


■/'*“!  / 


R c - 5c 
o 


dr, 


6T 


. /*  R 

1 / Is. 

~ ~ c I c 

o Jo  o 


dr . 


(3; 


In  expressions  (2)  and  (3),  the  6-operator  refers  to  differences 
between  different  realizations  of  the  acoustic  medium,  and  the  integrals  are 
taken  along  the  mean  ray  path.  The  phase  J associated  with  6T  is  <u6T  where 
u)  is  the  radian  acoustic  frequency,  therefore 


I 


_ u,  f R(6c 

cl  c 
oj  o o 


) dr . 


(4 


The  sound  velocity  perturbations  are  caused  by  vertical  shifting 
of  isovelocity  contours,  therefore  it  is  expected  that  6c  will  simply  be 
proportional  to  the  product  of  the  internal  wave  amplitude  t and  the  back- 
ground sound  velocity  gradient;  that  is 


r 3c 
6c  c o 

— « — 

c c 3z 
o o 


(5) 


Although  this  relation  is  conceptually  correct,  it  is  shown  by  Munk  (1974) 
that  the  gradient  3c  /3z  should  actually  be  replaced  by  the  gradient  of 


potent  ia 1 sound  velocity  3c^/3z  where  c^  does  not  include  pressure  effects 


or  the  effect  of  the  adiabatic  lapse  rate. 


J 
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The  corrected  version  of  (5)  can  be  substituted  in  (4)  giving 


i(R)=  -!  / 


( r ) ) dr 


(6) 


The  mean-square  phase  difference  between  two  spatially  separated 
receivers  is  given  by 


<j>2=  <[j(r2 ) - i(Ri ) ]2>, 


which,  by  (6),  become^ 


A2  - — 

V ~ c2 
O 


r Ri 

// 


9Cp(x,  ) 3Cp(x2 ) 


c 3z 
o 


c 3z 
o 


< t(x,)c(x2)  > dr,dr: 


3c  (x.  ) 3c  (x2  ) 

— E P 


c 3z  c 3z 

o o 


< c(x, )c(x2  ) > dr,dr2 


(7) 


R1R2 


_ 2 


3c  (x.)  3c  (x2) 

P P 


c 3z  c 3z 

o o 


< t(x, )c(x2 ) > dr ! dr j 


In  order  to  determine  typical  values  of  <J>2  , several  assumptions 


will  now  be  invoked,  namely,  straight  acoustic  rays,  constant  3c  /c  3z,  and 

P o 


stationary  < c(xi)c(x2)  > given  in  terms  of  the  Garrett-Munk  spectrum.  The 
assumption  of  straight  acoustic  rays  can  be  relaxed  somewhat.  Desaubies 
(1976b)  shows  that  ray  curvature  per  sc  is  mostly  unimportant  in  treating 


phase  perturbation;  however,  the  assumption  of  constant  Sc^/c^Sz  and 


uniform  <C2>  (in  depth)  are  important  and  are  not  trivially  removed.  Further 
comments  are  given  in  the  discussion. 


The  derivation  of  the  final  expression  for  $ 2 is  given  in 


Appendix  B;  the  result  is 


4>2  = <4tt 


A ( 90° ) 3Cp  r 

A . c 3z  Jrms 

acoustic  o 


( 8a  > 


T\ 


3c 

(“_)*(_& 
c 

o 
A 


. )2  - E b2(— ) A 2 ( 90° ) 
c 3 z it  n 

o 


R 


A ( 9 , ) ’ A (90°) 


, p , 0 j ,62! 


(8b) 


where  the  terms  E and  b come  from  the  internal  wave  spectrum,  R is  the 
median  source-receiver  range,  n is  the  Brunt-Vaisala  frequency  just  below 
the  thermocline,  n is  the  Brunt-Vaisala  frequency  at  the  transmission  depth, 
p is  n/n^,  n is  the  local  inertial  frequency,  A is  the  receiver  sepa- 
ration distance,  0t  is  the  angle  between  the  array  axis  and  the  median 
propagation  direction,  02  is  the  angle  between  the  array  axis  and  vertical, 
and  A(92)  is  the  relevant  internal  wave  length  scale,  which,  because  of 
vertical  anisotropy  in  the  Garrett-Munk  spectrum,  is  dependent  on  the  angle 
02.  The  function  T contains  the  entire  dependence  of  <j>2  on  receiver  sepa- 
ration, A,  receiver  orientation,  and  source-receiver  range,  R.  This 
dependence  is  illustrated  in  Figure  2 for  a horizontal  array  in  endfire 
orientation  (0,=O°,  02=9O°),  in  Figure  3 for  a horizontal  array  in  broadside 
orientation  (9]=90°,  02=9O°),  and  in  Figure  4 for  a vertical  array  in  broad- 
side orientation  (91=90°,  02=9O°).  In  Figure  2 there  is  no  range  dependence 

each  curve  refers  to  a different  value  of  p(=n(z)/n,  ) as  indicated.  In 

in 

Figures  3 and  4,  A ^ is  A(90°),  A^  is  A(0°),  each  group  of  curves  pertains 
to  a different  range  (log  R/A  is  given  for  each  group),  and  within  each 
group,  each  curve  pertains  to  a different  p (as  shown  on  right  side).  For 


LOG  1 0 (DELTR/LAMBDAH) 


Figure  2 


Horizontal  array,  endfire.  The  function  T is  used 
in  Eqs.  (8a)  and  (8b).  Seven  curves  are  shown;  each 
pertains  to  a different  value  of  u,  and  the  smallest 
and  largest  p-values  used  are  indicated  on  the  right 


LOG  1 0 (DELTA/LRMBOflH) 

Figure  i.  Horizontal  array,  broadside.  The  function  T is  used 
in  Eqs.  (8a)  and  (8b).  Values  of  log10(R/A[|)  are 
given  for  each  group  of  curves.  Within  each  group, 
seven  curves  are  shown,  each  for  a different  value 
of  ii,  and  the  smallest  and  largest  p-values  used  are 
indicated  on  the  right  for  one  of  the  groups.  The 
crosses  indicate  separations  at  which  the  acoustic 
phe sc  correlation  is  0.5  (not  to  be  confused  with 
acoustic  signal  correlation). 
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small  A/A(9),  T is  proportional  to  A2  for  the  horizontal  endfire  case,  and 


porpor t iona  1 to  A2  In  A for  the  other  two  cases.  For  very  large  A / A ( 9 J , T 


is  ultimately  proportional  to  Al n A for  the  horizontal  cases  and  proportional 


to  A for  the  vertical  case. 


The  internal  wave  length  scale  A(B)  is  given,  for  the  horizontal 


and  vertical  array  cases  respectively,  by 


x _ 1 con'll  - h n°  tan  1 ^n2-  T - 'V-  1 / m 2 ( 

H ' ; b 1 2n  In  p - (u2-  1 ) / 2 u 2 

in 


A = A (0° ) = b f ez/b 

V 3tt 


or,  since  u>>  1 usually, 


~ — 0 1 

H 8 n . ln|j  - \ ’ 
in 


b z/b 
X ~ — p 
V 9 e 


In  (8b)  the  parameter  E is  a dimensionless  "universal  constant";  it  repre- 


sents the  overall  internal  wave  spectral  level  and  it  has  a value  of 


5.3  x 10  . The  parameter  b represents  the  exponential  fall-off  rate  at 


large  depths  for  the  mean  Brunt-Vaisa la  frequency.  In  the  Garrett-Munk 


internal  wave  spectral  model,  the  Brunt-Vaisala  frequency  n(z)  is  taken  to 


n (z ) = n e 
o 


in  the  main  body  of  the  ocean  (i.e.  below  the  thermoc 1 ine ) . 


3- 


It  can  be  seen  from  Figures  2 to  4 that  the  rate  of  increase  of 
T with  A begins  to  reduce  near  A « A(0)  (for  fixed  R/A^).  This  occurs 
because  the  correlation  between  the  phase  parts  of  the  signals  at  the  two 
receivers,  which  is  almost  unity  for  A < A f © ) , begins  a more  rapid  decrease 
as  A increases,  and,  as  A becomes  larger  chan  A(0),  the  correlation  becomes 
significantly  less  than  inity.  The  phase  correlation  distance,  that  is,  the 
value  of  A at  which  the  correlation  of  the  phases  is  0.5,  is  not  reached, 
however,  until  A/A(e)  has  increased  to  the  value  indicated  by  the  cross 
(for  each  R/A^  curve). 

For  the  vertical  array  configuration,  T actually  decreases  as  A 
increases,  for  A/A^~  1.  More  specifically,  this  behaviour  exists  for  a 
particular  range  of  ray  angles,  measured  with  respect  to  the  horizontal, 
namely  10°  ~ I ray  angle|  £.85°,  0r  1 < (A/RKA^/A^  X.  100.  The  reduction 
occurs  because  is  significantly  greater  than  A^.  For  almost  horizontal 

propagation,  the  phase  coherence  length  is  controlled  by  the  horizontal 
internal  wave  coherence  length  A . As  A increases  with  R fixed,  the 
effective  coherence  length  tends  towards  A^,  i.e.,  it  reduces  significantly, 
with  the  result  that  significantly  more  of  the  phase  fluctuation  cancels 
itself  out  when  integrated  along  the  total  ray  path.  For  A >>  A(6),  T again 
increases,  but  this  is  because  in  this  regime  the  actual  propagation  range 
/R2+  \ A2  , significantly  increases  with  A,  along  each  constant  R/A  curve, 
and  it  more  than  compensates  for  any  reduction  tendencies. 

The  parameter  A ( 0 ) thus  serves  to  divide  the  separation  dis- 
tances into  two  regions:  one  region  includes  the  very  tightly  correlated 
separations  where  T is  small  but  increases  rapidly;  the  other  includes  the 
more  loosely  correlated  separation  where  T is  more  nearly  its  maximum  and 
increases  much  more  slowly  (not  including  very  large  A). 
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The  range  dependence  of  T can  be  summarized  as  follows:  for 
the  endfire  case  T is  range  independent;  for  the  horizontal  and  vertical 
broadside  cases,  T is  proportional  to  R In  R for  R and  4 ♦ 0,  it  is  range 
independent  for  A/X(6)-»  <*> , and  it  is  proportional  to  R between  these  two 
limiting  cases.  For  the  endfire  case  the  two  rays  are  identical  between 
the  source  and  the  nearest  receiver  and  thus  the  phase  difference  generated 
over  this  part  of  the  propagation  path  is  zero.  The  entire  phase  difference 
is  due  only  to  propagation  between  the  two  receivers  and  hence  is  indepen- 
dent of  the  source-receiver  distance.  For  the  other  two  cases  the  range- 
independence  at  very  large  A/X( 9)  occurs  because  A >>  R is  the  regime  and 
so  the  total  propagation  path  length  ^R2^-  is  virtually  independent  of  K; 

in  the  more  interesting  regime  of  small  A/X(6),  the  R In  R or  R-dependence 
is  due  to  the  detailed  form  of  the  Garrett-Munk  internal  wave  spectrum. 

For  A/X(9)«l,  Figures  2 to  4 show  that 


f 0.4 

4/)h 

for 

endfire  horizontal, 

/T  « j 

10.4, 

1 

^1R/ \ 

for 

broadside  horizontal, 

f 

, 0.4' 

Ja_  _R_ 

W \ \ 

for 

broadside  vertical. 

i 

1 A A 

v H 

By  using  (11)  for  this  range  of  A,  and  by  taking  the  square  root  of  (8b', 

it  can  be  seen  that  the  rms  phase  difference,  rms  <f> , is  proportional  to  the 

potential  sound  velocity  gradient  3c  / 3c,  the  acoustic  frequency  f (=w/2n), 

P a 

the  square  root  of  the  oceanographic  length  scale  b and  either  A or  * AR 

z / 2 b 

depending  on  the  configuration.  It  is  also  proportional  to  the  factor  e 
and  a term  dependent  only  on  p . This  is  summarized  as  follows:  for 
A/X(9)ftS  1 and  f /c  : X , the  acoustic  wavelength, 


1 
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rms  $ 


"t  endfire  horizontal. 

Ka 


2b(r~^)  VHtan~‘  /^zri  ~ /pI-  1/ui  r 


broadside  horizontal, 


H ’/*I 

" X * 
Av  a 


broadside  vertical. 


where  rms  6 is  given  in  degrees.  The  product  of  the  first  three  terms  in 

(12)  is  C p / ^E , where  C is  the  local  rms  internal  wave  height, 

rms  — r / rms  6 

c 3z / 

o' 

and  it  has  been  shown  elsewhere  (e.g.  see  Desaubies,  1976b  Eq.  11  or  MZ  Eq.  89) 

z/’b  3C 

that  although  C is  expected  to  increase  with  depth  as  eZ  “ , — is 
° rms  r r c03z 

expected  to  decrease  with  depth  as  e (at  least  for  the  profile  given  in 

10).  This  leads  to  the  conclusion  that  rms  4>  is  expected  to  decrease  with 
, . . , -1.5z/b 

depth  approximately  as  e . Thus,  for  straight  ray  propagation,  more 

scattering  by  internal  waves  should  occur  for  shallow  paths  than  deep  paths. 

Curved  rays  are  treated  for  a single  receiver  by  MZ  who  give  an 

exhaustive  treatment  of  scattering  in  convergence  zone  propagation  (excluding 

reflections).  They  find  that  most  of  the  scattering  occurs  near  the  apex 

of  the  ray  paths.  The  reason  again  is  that  C3c  /3z  follows  a similar 

P 

exponential  depth-dependence  as  the  Brunt-Vaisala  frequency  (Eq.  10)  and  so 
the  internal  waves  produce  more  scattering  (possibly  orders  of  magnitude  more) 
near  the  thermocline  than  in  abyssal  regions. 

Desaubies  (1976a)  has  provided  an  interesting  and  very  attrac- 
tive variant  of  the  Garrett-Munk  spectrum.  He  shows  that  the  three  parameters 

E,  N and  b can  be  reduced  to  two  parameters  r and  t (using  his  notation) 

o 


r = E n b2 
o 


-16- 

t = 3/  ( 2bn  ) . (141 

o 

With  these  definitions,  the  dependence  of  n2  on  z need  not  be  exponential 

but  may  remain  unspecified  and  completely  general.  The  scaling  parameters 

in  the  spectrum  cease  to  be  E,  n , n.  and  b and  become  r,  t (neither  of 

o in 

which  is  dimensionless),  n.  and  the  local  frequency  n(z).  He  further  shows 

in 

that  the  existing  data  strongly  suggest  that  r and  t are  universal  constants 

-4 

with  values  of  300  m2cph  and  3.10  (cycles  per  meter)  /cph  respectively. 

This  is  equivalent  to  fixing  the  values  of  n^  and  b,  by  equations  (13)  and 

z / b 

(141  to  be  universally  4.42  cph  and  1132  m respectively,  and  replacing  e 
by  4.42/n(z)  where  n(z)  is  left  unspecified.  In  Appendix  A,  Equations  (8b), 
(9)  and  (12)  are  rewritten  using  the  r,  t scaling. 

TYPICAL  OCEANIC  PARAMETERS 

Numerical  estimates  will  now  be  obtained  for  four  different 
oceanic  regions,  see  Figure  5;  the  Gulf  of  Alaska  (North  Pacific),  the 
Newfoundland  Basin  (North  Central  Atlantic),  Robeson  Channel  in  April  (82°N, 
60°W  in  the  Eastern  Arctic),  and  the  Canada  Basin  in  April  (78°N,  130°W  in 
the  Western  Arctic).  In  each  region  estimates  are  obtained  for  an  acoustic 
frequency  f of  150  Hz,  a range  R of  50  km  and  a hydrophone  spacing  of 
either  1 km  for  the  horizontal  separations  or  50  m for  the  vertical  sepa- 
ration . 

Internal  wave  measurements  in  the  temperate  oceans  indicate  that 
E = 5.3  10  and,  in  the  absence  of  Arctic  measurements,  this  value  will  be 
used  for  the  Arctic  case  as  well.  Table  1 contains  all  the  numerical  values 
that  have  been  used.  The  Gulf  of  Alaska  and  Newfoundland  Basin  data  were 
calculated  from  temperature  and  salinity  profiles  published  by  Hashimoto 
(1968,  his  regions  D and  P respectively).  At  the  chosen  depth  (1000  m)  the 


Figure  5.  The  solid  regions  indicate  the  locations  of  the  four 
numerical  examples:  Gulf  of  Alaska  (North  Pacific), 
Newfoundland  Basin  (North  Central  Atlantic),  Robeson 
Channel  (82°N,  60°W),  and  Canada  Basin  (78°N,  130°W). 


r 
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profiles  are  independent  of  season.  The  Robeson  Channel  data  were  calcula- 
ted from  temperature,  salinity  and  density  profiles  published  by  Sadler 
(1976,  his  Figure  23)  and  the  Canada  Basin  data  were  taken  from  Herlinveaux 
(1463,  Station  1,  Camp  1).  All  n(z)  profiles  were  calculated  using  the 
formula  n2(z)  = (g/;)(3c/3z)  - g2/c2(z),  and  in  situ  density  values  were 
calculated  using  the  equation  of  state  given  by  Chen  and  Millero  (1976). 

The  values  of  n^  and  b given  in  Table  1 were  obtained  by  fitting  (10)  to 

flu*  relevant  profile;  values  of  n in  Table  I were  then  calculated  from  (10) 

o 

using  n , b and  the  relevant  depth,  z.  The  sound  velocity  profiles  were 
calculated  using  Wilson's  equation;  the  potential  sound  velocity  profiles 
were  calculated  by  first  determining  the  local  adiabatic  sound  velocity 
gradient  from  the  local  adiabatic  lapse  rate  (Fofonoff,  1962)  along  with 
the  local  salinity  gradient  and  0.  17<>  of  the  local  pressure  gradient  (this 
very  small  correction  more  properly  treats  the  internal  wave  pressure  field;, 
then  integrating  the  adiabatic  profile  from  the  surface  down  and  finally 
subtracting  the  resulting  values  from  the  total  sound  velocity.  Gradients 
in  the  potential  sound  velocity  were  obtained  by  numerically  differentiating 
the  final  profile  with  respect  to  depth. 

The  rms  phase  fluctuations  and  the  internal  wave  length  scales 
calculated*  from  these  data  are  given  in  Table  II.  The  most  striking  result 
is  the  comparatively  large  phase  instability  estimated  for  the  Newfoundland 
Basin.  It  is  a direct  result  of  the  relatively  large  potential  sound 
gradient  existing  in  that  region  which  in  turn  may  be  due  to  its  more 
southerly  location.  The  other  estimates  are  more  uniform,  with  the  Arctic 
Values  equivalent  to  Table  II,  but  using  Desaubies'  modified  scaling 


1 


! 


are  given  in  Appendix  A. 
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TABLE  1 

OCEANOGRAPHIC  DATA  USED  IN  THE  NUMERICAL  EXAMPLES 


Region  and  Depth 

b(m) 

3c  i 

rfe  <■  > 

o 

c (m/s) 

0 

n (cph) 
o 

n(cph ) 

n , (cph ) 
in 

Gulf  of  Alaska 

(1000  m) 

1900 

-3.51  10~6 

1477 

2 

1.18 

i 

15.7 

Newfoundland  Basin 

(1000  m) 

2 200 

-41.5  10-6 

1501 

. 

3 

1.90 

1 

17.0 

Robeson  Channel 

(Apri 1 , 300  m) 

86 

-6 

+3.95  10 

i 

1453 

30 

0.90 

1 

12.0 

Canada  Basin 

(April , 1000  m) 

1000 

-1.53  10-6 

1466 

2 

0.74 

1 

12.0 

TABLE  II 

RESULTS  FOR  A RANGE  OF  50  KM,  AN  ACOUSTIC  FREQUENCY  OF  150  HZ  AND  A 
RECEIVER  SEPARATION  OF  EITHER  1 KM  (HORIZONTAL  CASES)  OR  50M  (VERTICAL  CASE) 


Region  and  Depth 


Endf ire 
Hor  izonta 1 


Gulf  of  Alaska  (1000  m) 


Newfoundland  Basin  (1000  m) 


Robeson  Channel  (April  300  m) 


Canada  Basin  (April  1000  m) 


Horizontal  Broadside, 
rms  <t>  = < [«(R2)-K(Ri)]J>^ 


Broadside 
|Hor  izonta  1 


1.5 


20 


0.31 


0.40 


4.6 


53 


1.2 


1.5 


Broadside1 

Vertical 


3.9 


65' 


0.71 


0.86 


3.0  km 


4.7  km 


1.9  km 


340  m 


370  m 


300  m 
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being  more  stable  and  having  somewhat  shorter  internal  wave  length  scales. 
The  increased  stability  in  the  Arctic  is  due  to  weaker  internal  wave  acti- 
vity, as  predicted  by  the  Garrett-Munk  model;  the  gradients  of  c are,  in 
fact,  almost  the  same  as  the  Gulf  of  Alaska  value. 

Under  the  assumption  that  the  internal  wave  statistics  are 
Gaussian,  the  acoustic  signal  coherence  between  the  two  receivers  drops  to 
0.5  when  rms  <f  reaches  67.4°  (keeping  acoustic  amplitudes  fixedj.  Figure  6 
illustrates  the  dependence  of  receiver  separation  on  horizontal  median— ray 
range  for  an  acoustic  signal  coherence  of  0.5  for  each  of  the  four  oceanic 
regions.  Only  horizontal,  broadside  separation  is  shown,  and,  again  the 
frequency  is  150  Hz.  In  the  N.  Atlantic,  short  to  moderate  ranges  produce 
a coherence  of  0.5  (e.g.  50  km  range,  2.3  km  separation),  in  the  N.  Pacific 
the  ranges  and  separation  are  considerably  larger  (e.g.  60  km  range  and 
30  km  separation),  and  in  the  Arctic  the  range  and  separation  are  so  large 
that  they  are  probably  outside  the  framework  of  this  analysis. 

DISCUSSION  AND  CONCLUSIONS 

Internal  waves  cause  phase  decorrelations  between  signals 
transmitted  by  a single  source  and  received  simultaneously  by  two  separated 
hydrophones.  The  rms  phase  difference  is  a rapidly  increasing  function  of 
the  separation  distance  if  the  separation  is  less  than  the  internal  wave 
length  scale,  2 to  3 km  for  a horizontal  array  or  350  m for  a vertical 
array*.  The  rms  phase  difference  increases  more  slowly  for  greater  sepa- 
ration distances.  At  150  Hz  and  50  km  range,  the  rms  phase  difference  is 
expected  to  be  quite  small  for  both  the  Arctic  and  the  Gulf  of  Alaska,  but 


it  is  expected  to  be  much  larger  in  the  Newfoundland  Basin,  mainly  because 
* Using  Desaubies  scaling,  these  numbers  become  3 to  4 km  for  a horizontal 
array  and  300  to  700  m for  a vertical  array. 


RECEIVER 


the  product  of  the  potential  sound  velocity  gradient  and  the  rms  internal 
wave  height  is  much  larger  there  than  in  the  other  regions.  If  the 
frequency  is  increased,  the  rms  phase  differences  will  increase  propor- 
tionately; if  the  range  is  increased,  the  rms  phase  differences  will  again 
increase  but  only  as  the  square  root  of  the  range  (approximately).  In 
either  case,  the  internal  wave  length  scales  will  increase.  For  conver- 
gence zone  propagation,  it  has  been  shown  elsewhere  that  most  of  the  phase 
difference  is  produced  at  the  shallowest  parts  of  the  ray  paths. 

In  order  to  provide  an  appreciation  for  some  of  the  effects  of 
these  phase  perturbations,  three  very  simple  calculations  have  been  made, 
each  using  a horizontal  broadside  array  in  the  Newfoundland  Basin.  The 
three  calculation  are:  array  gain  loss,  beam  steering  error,  and  range 
error  using  wavefront  curvature  estimates.  For  an  array  with  M hydrophones 
operating  in  a non-perturbed  medium,  the  array  gain  in  dB  is  10  log(M)  if 
the  signal  is  perfectly  coherent  and  the  background  noise  is  incoherent 
between  hydrophones.  The  gain  is  significantly  reduced  in  the  presence 
of  phase  perturbations  if  the  array  is  much  longer  than  the  acoustic  cohe- 
rence length  as  given  in  Figure  6.  In  the  present  case,  internal  waves 
reduce  the  gain  by  3 dB  for  an  array  length  of  approximately  5 km*.  For 
larger  arrays,  the  reduction  is  more  dramatic,  in  fact,  the  gain  ultimately 
becomes  independent  of  array  length.  It  should  be  noted  that  this  loss 
is  frequency  and  range  dependent;  the  loss  will  increase  as  the  range 
increases  and  decrease  as  the  acoustic  frequency  decreases.  It  should  also 
be  noted  that  'approximately  5 km'  is  only  coincidentally  the  same  as  the 
* For  a horizontal  endfire  array  the  necessary  array  length  for  3 dB  loss 
is  approximately  12  km;  for  a vertical  broadside  array,  it  is  approximately 


2 50  m. 
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the  i n r o rna  1 wavf  coherence  length,  4.7  km.  The  beam  st  eering  error  is 
very  small.  Even  if  only  two  hydrophones  are  used,  the  rms  steering  error 
is  less  than  0.1  for  a hydrophone  separation  greater  than  700  m,  and  the 
error  decreases  as  the  separation  increases.  It  is  expec  ted  that  less 
error  would  result  if  a statistical  estimation  method  were  employed  using 
more  than  two  hydrophones.  This  error  is  frequency  independent  and  it 
increases  approximately  as  the  square  root  of  the  range.  The  range  esti- 
mation error  using  wavefront  curvature  is  also  frequency  independent  and 
decreases  as  the  hydrophone  separation  increases;  however,  it  increases  as 
approximately  the  3/2's  power  of  the  range.  At  least  three  hydrophones  are 
necessary  to  form  the  estimate,  and,  for  the  present  case,  with  equal 
hydrophone  spacings  of  0.5  km,  the  rms  range  error  using  only  3 hydrophones 
is  157,  of  the  range.  If  more  than  three  hydrophones  are  used  with  a 
statistical  fitting  technique,  this  error  will  probably  decrease. 

The  model  that  has  been  presented  here  is  highly  simplified  and 

rather  loosely  applied,  partly  because  it  assumes  that  the  Garrett-Munk  (or 

. 

Desaubies)  internal  wave  spectrum  applies  equally  well  everywhere.  There  is 

evidence  that  it  applies  in  the  temperate  oceans  (Desaubies,  1976a);  as  yet 

no  evidence  has  been  published  concerning  its  applicability  in  the  Arctic. 

The  other  important  simplification  is  that  the  potential  sound  velocity 

gradient  and  the  rms  internal  wave  height  are  depth-independent.  More 

specifically,  the  assumption  is  that  the  product  (3c  /c  3z)2A2  t2  is 

r p o H rms 

constant  along  an  acoustic  ray,  and,  for  medium  to  long  ranges,  this  is 
equivalent  to  the  previous  statement.  The  variation  of  this  product  is 
largely  controlled  by  the  term  ( 3c^/3z ) 2 /n (z ) and  this  in  turn  is  controlled 
by  the  local  potential  temperature  gradient  and  the  local  salinity  gradient. 


Munk  (1974,  his  equation  ( 1 0 1 ) shows  that  3c  / 3z  and  n2(z)  are  proportional; 
he  also  shows  that  the  "constant"  of  proportionality  is  a function  of  the 
Turner  number,  that  is,  a number  which  is  proportional  to  the  ratio  of  the 
local  salinity  and  potential  temperature  gradients.  The  term  (3c  / 3z j! /n(z  ' 
can  therefore  be  expected  to  vary  with  depth  and  with  geographical  location. 
As  shown  earlier  for  the  Garrett-Munk  profile,  it  varies  approximately 
with  depth  as  e and  from  this  it  can  be  seen  that  it  reduces  by  a 

factor  of  20  if  the  ray  depth  increases  by  an  amount  equal  to  b.  This 
reduction  factor  will  be  at  least  partly  offset  by  the  fact  that  the  rays 
must  necessarily  be  non-horizontal  over  part  of  the  path,  probably  that 
part  nearest  the  surface,  and  thus  for  the  "horizontal  configuration"  cases, 
scattering  will  be  partly  due  to  the  vertical  structure  of  the  internal 
wave  field  at  depths  when  the  vertical  structure  is  strongest.  This  struc- 
ture possesses  shorter  coherence  lengths  and  can  produce  much  larger  phase 
differences  than  the  horizontal  structure  (c.f.  Table  II  and  note  that  the 
vertical  receiver  separation  is  only  50  m whereas  the  horizontal  receiver 
separation  is  i000  m).  The  overall  effect,  then,  of  this  simplification 
is  in  fact  unknown,  but  it  is  expected  that  it  will  result  ir  estimates  of 
<f2  that  are  too  large. 

Removal  of  other  simplifications  in  the  work  presented  here  (e.g. 
only  three  array  configurations)  or  extensions  to  other  geographical  regions 
could  be  done  straightforwardly. 
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APPENDIX  A 


DESAUBIES ' SCALING  FOR  THE  GARRETT-MUNK  SPECTRUM 
Using  Desaubies'  scaling,  Eq.  (8b)  becomes 


2 9 c „ 

a)  . , p , , 2r 


*2  = (7L)  ( 90°) } 2 { tan-'  v/i72"-  F -'r^~l/V2}T,  (A1 ) 

c c d z Tin 
o o 


or,  with  u = °°  and  r = 300  m2cph, 


rms  <J>  = 6235  ( — •—)  - — l — degrees 
c dz  X n 
o a 


where  T = T at  p = n is  expressed  in  cph,  and  SI  units  are  used  else- 

ao 

where.  Also,  from  (9a)  and  (9b) 


*H  8tn. 

in 


1 1 

v 2irtn . p 
in 


1 tan  ^p 2 — 1 -^p2-  1/p2 


In  p - (p2-  1 ) / 2 m 1 


With  p >>  1 and  t = 3 10  (cycles  per  metre)/cph, 


A„  = 654.5/(n.  In  p), 
H in 


A = 530. 5/n 

V 


Table  Al  shows  the  results  equivalent  to  Table  II  but  calculated 
using  the  local  measured  value  of  n(z)  and  using  (Al,  A3  and  A4 ) . The  main 
difference  from  the  Table  II  values  occurs  for  the  Arctic  case;  Desaubies' 
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scaling  predicts  more  phase  fluctuation  and  longer  internal  wave  length 
scales  than  the  Garrett-Munk  scaling.  In  fact,  according  to  Table  Al,  the 
Arctic  and  Gulf  of  Alaska  display  essentially  the  same  degree  of  instability. 
The  Newfoundland  Basin  is  again  predicted  to  be  the  most  unstable  but  this 
time  with  the  shortest  vertical  internal  wave  length  scale  instead  of  the 
longest . 


A-3 

• • 

TABLE  a 1 

DESAI'BIES'  MODEL-  Phase  difference  for  a range  of  50  Km,  an  acoustic  fre- 


quency of  150  Hz  and  a receiver  separation  of  either  1 km  (horizontal  case) 
or  50  m (vertical  case). 


Region  and  Depth 

n (cph  ) 

Horizontal  Broadside  . 
rms  if  =<[l>(R2)-4>(Ri  )]2>^ 

XH 

\ 

Endf ire 

Hor izonta 1 

Broadside 
Hor izonta 1 

Broadside 

Vertical 

V 

Gulf  of  Alaska 
(1000  m) 

1.20 

1.3° 

3.9° 

3.5° 

4.0  km 

450  m 

Newfoundland  Basin 
(1000  m) 

2.55 

11° 

31° 

39° 

3.6  km 

280  m 

Robeson  Channel 
(April , 300  m) 

1.97 

1.1° 

m 

2.0° 

3.7  km 

590  m 

Canada  Basin 
(April,  1000  m) 

.66 

.74° 

2.2° 

1.1° 

4.0  km 

720  m 

APPENDIX  B 


MATHEMATICAL  DETAILS  LEADING  TO  THE  FINAL  EXPRESSION  FOR  $2. 


In  F.q.  (7)  on  page  6 the  term  <c(x1)c(x2)>  is  the  correlation 
function,  R^,  of  the  internal  wave  spectrum.  This  in  turn  is  the  Fourier 
transform  of  the  power  spectrum  in  wave  number  space,  and  it  can  be  shown 
that  the  slant  range  correlation  is  given  by 


R (H,V) 


r n(z)  /•  » 22. 

= / doi  / da  Z2  E(a,w)  cos[ 2waV(-^-j — j— ) 2]  J (2waH) 

I I a.  - n:  -1  o 

Jn  .Jo  in 


w2-  n2 
in 


where  Z2  = b2n  — — j — , E(a,u>)  = - EA(X)  113 , A(X)  = (2/tt)(1  +X2)" 

o nuj  7i  ✓ 9 « \ ^ 7 

u)(w2-  n2  )o 
in  * 


X = a/a^,  = (3/2nQb)  'u2-  n2^,  H = horizontal  component  of  xr  x2  , V is 

the  vertical  component  of  x,-  x2 , and  a is  expressed  in  cycles  per  metre 
(see  Desaubies  1976a,  Eq.  10,  and  note  that  R (H,V)  as  given  in  (Bl)  is  the 

CD 

same  as  hisf  dw  MSC(a))  MS(u>)).  By  substituting  this  expression  in  (7)  and 
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taking  the  gradient  — outside  the  integral  signs  (because  they  are  assumed 
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constant  over  the  propagation  path),  the  mean-square  phase  difference  becomes 
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In  this  expression,  V and  H are  function  of  r,  and  r 2 , but  their  functional 
forms  are  not  necessarily  the  same  for  each  of  the  three  double  integrals 
in  r,  , r2 . For  the  horizontal  endfire  configuration,  they  are  all  the  same: 

V = 0,  H = | r 2 — r 2 I , with  Kj  = R - 4/2  and  R2  = R + 4/2;  however,  for  the 

horizontal  broadside  configuration,  V = 0,  H = r,  - r2  for  the  first  two 
integrals  and  H = ^(r,-  r2  ) 2 + (r,+  r2)2  4 2 / 2 } / { 1 + 4 2 / 4R 2 } , with 

, = R2  = ^R2  + 42/4,  where  4 is  the  hydrophone  separation  distance.  For 
the  vertical  broadside  case,  again  the  first  two  are  the  same, 

V = \ |r2+  r,|  4/  >rR*  + 42/4,  H = |r2-  r ! | R/  ^R2  + 42  /4 , and  the  third  is 

V = (4/R)  (r2±  k /(r2-  r,)!  + 2(r?  - rj)  42/4R2+  (r2  + r,)2  4-/16R4) 

X( 1 +42/4R2)"3/2,  H = (r242/2R2  ± /(r2-  r,)2+2(r?-  r 2 )4 2 /4R2+(r  2 ♦ r ~ , 2 4 V 1 6R4  . 

X(  1 +42/4R2)  ^2 , with  R(  = R2  = ^R2  + 4 2 /4 . The  general  form  given  in  (B2) 

can  be  further  simplified  by  replacing  J (z)  with  2/ttP  ^ cos  (zs  in6  ;df . 
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The  two  cosine  terms  have  arguments  that  are  linear  in  u and  they  can  be 
combined.  The  integration  in  u can  be  carried  out  with  the  help  of  the 
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formulaj  du  cos(ua)/(l  + u2  ) = tt/2  e 
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dr , dr 2 


n -3irn.  
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-|V  ^y2-  x2  - H "x2-  1 sin9| 


In  all  cases,  the  integrals  in  r,  and  r2  can  be  carried  out,  but,  because 
of  algebraic  complexities,  only  the  horizontal  endfire  case  will  be  shown 
here.  For  this  case,  V = 0 and  dr,  dr2  can  be  converted  to  dr  dH  with 
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appropriate  sets  of  limits.  The  integrals  can  then  be  done  by  parts  and 
the  final  result  is 


n/2  p 
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— r—  A^x2 - 1 sin9  3nn  . 

n b in  r~z 

o + r—  A'x2-  1 sine  - 1 
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The  combination  of  parameters  n^b/Bun^  is  proportienal  to  (see  Eq.  9a  i, 
and,  in  the  limit  A * o , the  double  integral  becomes  — A2{tan  * ''p2-  F 

TT 

- 2 - \ / v2 ) . Introducing  these  terms  explicitly  into  (B4)  produces 


t7  (rf^)2  * Eb*("ir)  ah  {tan_1/^- 1 - -v-  i/p2)t, 
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? t C U ->  T~  P(u) *^x 2 — 1 Sine 

- I del  dx  1 i e H + P(  u ) ^ x2  - T sin6 

'J0  J t P2(u)x’(x2-  1 )*s in2  6 H 


tan  *m2-  1 -'p!-  1 / m 2 


P ( M ) = - (tan  ' ^ m 2 - 1 -^p2-  l/p2)/(lnp  - (p2-  1 ) / 2 p 2 ) (B7) 


The  transformation  just  introduced,  and  the  resulting  equation  (B5),  are 
also  valid  for  the  general  case.  For  the  endfire  configuration,  it  can  be 


seen  that  T is  independent  of  range,  and,  as  A-*o,  T*fc(A/Xu)2,  independent 

of  p.  (As  A-*’®,  T-»  (^7  In  p)  In  —•) 
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In  general,  the  rms  internal  wave  height  c can  be  obtained 

rms 

from  (Bl)  at  H = V = 0: 


tions,  and  the  result  is  (8a).  The  general  form  of  T can  be  obtained  by 
comparing  (B5)  and  (B3): 


The  functional  description  of  T in  (8b)  is  equivalent  to  ( B8 ) . 


